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Abstract 

There currently exists no algebraic algorithm for computing twisted 
conjugacy classes in free groups. We propose a new technique for deciding 
twisted conjugacy relations using nilpotent quotients. Our technique is a 
generalization of the common abelianization method, but admits signifi- 
cantly greater rates of success. We present experimental results demon- 
strating the efficacy of the technique, and detail how it can be applied in 
the related settings of surface groups and doubly twisted conjugacy. 

1 Introduction 

Given two elements x and y in a group G and an endomorphism ip : G — > G, we 
say that x and y are twisted conjugate if there is some z such that 

x = (p(z)yz^ 1 . 

Twisted conjugacy is a generalization of the ordinary conjugacy relation in 
groups, and the computation of twisted conjugacy classes is a problem of con- 
siderable difficulty for many groups G. 

Computing twisted conjugacy classes (also called "Reidmeister classes") is 
of interest in various algebraic contexts. Our own approach will be motivated by 
Nielsen fixed point theory, though other motivations exist (see [Bogopols ki et ai, 2006| , 
in which the twisted conjugacy problem in free groups arises naturally in the 
context of the ordinary conjugacy problem in certain other groups). 

Extant algebraic techniques for computing twisted conjugacy classes are ad 
hoc in nature, and the goal of this paper is to present a new technique which is 
more generally applicable (though still not in general algorithmic), along with 
experimental results demonstrating its success. 
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Our technique is an extension of the common abelianization technique. View- 
ing the abelianization as the first nilpotent quotient, we show that twisted con- 
jugacy classes can be distinguished with increasing success rates when projected 
into nilpotent quotients of increasing nilpotency class. These projections, and 
the computations necessary to compute twisted conjugacy in nilpotent groups, 
can be fairly intensive, and as such we have implemented our technique in the 
Magma programming language. 

In Section [2] we give the motivation for the twisted conjugacy problem from 
Nielsen theory, Section [3] is an outline of our technique, Section 2] gives some 
examples, Section [5] presents our experimental results, and Sections [5] and [7] 
show how the technique can be adapted into two natural generalizations of the 
main problem. 

The bulk of this work was completed as part of the author's doctoral disser- 
tation, advised by Robert F. Brown. The content in this paper was also advised 
and suggested by Peter Wong. The author wishes to acknowledge their help 
and support, as well as helpful input from Seungwon Kim, Armando Martino, 
and Evelyn Hart. 

2 Nielsen fixed point theory 

Our principal motivation for studying the twisted conjugacy problem is Nielsen 
fixed point theory (standard references are |Jiang, 1983| and |Kiang, 1980| ). 
Given a map / : X — > X of a space with universal covering space X and 
projection map g : X — > X^the fixed points of / are partitioned into classes of 
the form p(Fix(/)), where / is a lift of /. If a single lift / is fixed once and for 
all, each fixed point class can be expressed as p(Fixa -1 /) for some a <E iri(X). 
A central problem in Nielsen fixed point theory is to determine the number of 
"essential" fixed point classes of a mapping. This number gives a lower bound 
for the minimal number of fixed points for maps in the homotopy class of /, and 
in the case when X is a manifold of dimension at least 3, this "Nielsen number" 
is in fact equal to the minimal number of fixed points. 

A fundamental question when counting fixed point classes is the following: 
given two elements a, f3 £ tt\(X), when is p(Fix(a -1 /)) = p(Fix(/3 -1 /))? The 
question can be answered with an elementary argument in covering space theory. 
Two such fixed point classes are equal if and only if there is some 7 £ tti(X) 
with 

a = </?(7)/37 -1 , 

where tp : tti(X) — > tti(X) is the map induced by / on the fundamental group. 
If the above holds, we say that the elements a and [3 are twisted conjugate (if 
if is the identity map, this is the ordinary conjugacy relation). The twisted 
conjugacy classes of / are also called Reidemeister classes, and the set of such 
classes is denoted Tl(<p). 

In the case where / is a selfmap of a compact surface with boundary, the 
fundamental group t^i(X) will be a free group, say ^i(X) — (gi, . .. ,gk)- No 



2 



general algorithm is known for computing the Nielsen number of such a map (the 
case k = 1 is classically known, and an algorithm for the case k = 2 has recently 
been developed by Yi and Kim |Yi and Kim, 2008] , extending work by Joyce 
Wagner Wagner, 1999 ). Fadell and Husseini, in |Fadell and Husseini, 1983] , 
proved: 

Theorem 2.1 (Fadell, Husseini, 1983). Let tp be the map induced by f on 
tti(X). Then the Nielsen number of f is the number of terms with nonzero 
coefficient in 

( n p, 

where d is the Fox Calculus operator (see \Crowell and Fox, 1963^ ), and p : 
Z7Ti(Y) — > UZ(ip) is the linear extension of the projection into twisted conjugacy 
classes. 

The theorem above thus reduces the problem of computing the Nielsen 
number to the application of the projection p, which requires some algorithm 
for distinguishing twisted conjugacy classes. No such algorithm exists in the 
literature, though Bogopolski, Martino, Maslakova, and Ventura give an al- 
gorithm in |Bogopolski et al., 2006| for the case where ip is an isomorphism. 
Their algorithm involves the traintracks machinery of Bestvina and Handel 
|Bestvina~a nd Handel, 1992 . This paper explores a purely algebraic technique 
which can be implemented by existing computer algebra systems, and in some 
cases can be done by hand. 



3 Abelian and nilpotent quotients 

Throughout this and the next two sections, let G be a finitely generated free 
group, and let ip : G — ► G be an endomorphism. For an element g G G, let 
[g] denote the twisted conjugacy class of g. Our goal is to confirm or deny the 
equality [g] — [h] for two group elements g and h. 

Existing algebraic techniques for deciding equality of such classes are sur- 
veyed in |Hart, 2005] . One such technique is the algorithm of Wagner Wagner, 1999 , 
which is applicable if the mapping tp satisfies the combinatorial "remnant" con- 
dition (this condition is satisfied with probability approaching 1 as the word 
lengths of the images of generators increase). For general mappings (without 
remnant), distinguishing classes is typically done by projecting into the abelian- 
ization G and solving the twisted conjugacy relation there. 

Abelianization is often sucessful at showing that two classes are distinct, but 
cannot be used to show that two classes are equal. Projecting into nilpotent 
quotients is a natural generalization of the abelianization technique, and we 
show how it can also furnish a technique for equating classes. 

We will use the commutator notation [x,y] = xyx~ x y~ x . Let 7n(G) be 
the terms of the lower central series, and let G n = G/^ n {G). Each of the 
groups G n are nilpotent of class n. The abelianization is of course G = G\ = 
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G/ji(G). In general, we will use a bar to indicate projection of elements into 
the abelianization, and a hat to indicate projection into G n , with the value of 
n to be understood by context. 

Computation in G is made easy by commutativity. For n > 1, the groups 
G n are not commutative, but powerful commutation rules make computation 
possible. The following easily verifiable commutator identities hold in any group: 

yx = [y,x]xy, 
[y,x] = [x,y]~ l , 
[xy,z] = [x, [y,z]][y,z][x,z] 

These rules have a nicer form in a class 2 nilpotent group, where all com- 
mutators will freely commute: 

Proposition 3.1. If G is a class 2 nilpotent group, then, for any x,y,z G G, 
we have 

yx = xy[x,y]' 1 , (1) 

[y,x] = [x,y]-\ (2) 

[xy,z] = [x,z][y,z]. (3) 

In a class 2 nilpotent group, we may use |T]) to exchange the order of any non- 
commutator elements. Using (J3]) we may write any commutator as a product of 
commutators of generators. Having reduced all commutators to commutators of 
generators, we may use to ensure that the generators appear in a prescribed 
order. Viewed as a set of word rewriting rules, Proposition 13.11 suggests that 
there will be some sort of normal form in nilpotent groups which can be used 
to compare words. 

The desired normal form is provided by a theorem of P. Hall. Before stating 
the theorem, we give some terminology and notation, following |Hall, 1957] . 

For a free group G, consider the elements which can be formed by taking the 
closure of the generator set under the commutator operation. Of these elements, 
the generators are referred to as weight 1 commutators, and the weight of any 
non-generator element is defined to be the sum of the weights of the elements 
of which it is a commutator. 

Hall showed that words in G n can be given in a normal form consisting of a 
product of certain basic commutators of weight n or less given in some proscribed 
order. The construction of the basic commutators is somewhat involved, and we 
refer to [Hall, 1957 and |Magnus et at, 19 76 for the details. For the purpose 
of the examples in this paper, it is sufficient to know that in a group of rank 
2, say G = (a,b): the basic weight 1 commutators are a and b, and the only 
basic weight 2 commutator is [a, b] (the commutator [b, a] is not basic, since it is 
expressable as [a,b}~ 1 ). We also refer to Theorem 5.11 of Magnus et ai, 1976| 



which gives a combinatorial formula due to Witt for C n , the number of basic 
weight n commutators: 

n 

d\n 
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where [i is the Mobius function, and k is the number of generators of G. 

Theorem 3.2 (P. Hall, 1957). For any x G G, we can write the projection 
x G G n as 

* = 11^' 

i 

where {cj} is the sequence of basic commutators of weight less than or equal to 
n. 

This form for x is unique up to the ordering of the weight n basic commu- 
tators, and we call this the Hall normal form for x. 

Since <^(7„(G)) C 7„(G), there is a well defined quotient mapping <p : G n — > 
G n . Thus it makes sense to ask, for h,gdG, whether or not [h] = \g\ in G n , 
that is, whether or not there is some z G G n with 

h = (p(z)gz~ l . 

If [h] 7^ \jg] in G n , then we know that [h] ^ [g] in G. 

4 Some examples 

We begin with a sample computation by hand, showing how the Hall normal 
form can be used to solve twisted conjugacy relations. 

Example 4.1. We will compute the Nielsen number of the map on a surface 
with fundamental group G — (a, b) which induces the homomorphism: 

a i— > ab 

^ b 

Theorem 12.11 gives 

RT(<p) = p(-l-b) 

and thus we need only decide whether or not [1] = [b]. First we attempt to 
equate these classes in the abelianization G. Writing elements additively, an 
element z G G is of the form z — na + mb, and we wish to solve 

= Cp{z) + b — z. 

We compute (p{z) = n(a + b) + m(4a + 26) and — z = —na — mb, and the above 
equation becomes 

—6 = 4ma + (n + m)b, 
and we can solve for n and m to find that z — —a is a solution. 

Remark 4.2. Given that 1 and b are twisted conjugate in the abelianization by 
the element —a, we might hope that 1 and b are twisted conjugate in the group 
G by the element a -1 . This is not the case, however, as tp(a~ 1 )ba = b~ 1 a~ 1 ba. 
The possibility remains, however, that these elements are twisted conjugate by 
some more complicated word which abelianizes to —a. 
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Having failed to decide the twisted conjugacy after a check in the abelian- 
ization, we proceed to the class 2 nilpotent quotient G2. Any z S G2 is of the 
form z — a n b m \a, b] k . We wish to solve 

1 = 

We already know by the above calculation, however, that any such element z 
must abelianize to a £ G in order to satisfy the twisted conjugacy equation. 
Thus we may assume that z — a -1 [a, b] k . Now we compute 

z^ 1 = [a,S] _fc a = a[a, b]~ k , 

and 

0(z) =6- 1 a- 1 [a6,6 2 ^] fc = a- 1 &- 1 [a-S&- 1 ][a,6] 2fe [6,a] 4fe = a" 1 ^ 1 ^, 6]- 2fe+1 , 

where we have used the rules of Proposition 13.11 together with the identity 
which follows from setting x — y in identity ([3]). 
We are now ready to test the twisted conjugacy equation above. The right 
hand side is: 

^(z)Sz- 1 ^a- 1 b- 1 [a,^~ 2k ~ 8 ba[a,b]- k = [a,b]~ 3k ~ 8 . 

Setting this equal to 1 requires that —3k — 8 = 0, which is impossible for fceZ. 
Thus there can be no such z 6 G2, and so [1] ^ [b] and the Nielsen number is 2. 

The example above involved a computation first in G and then in G2 ■ In each 
step, there are essentially two types of computational operations involved. The 
first is the term rewriting to obtain the Hall normal form, which was fairly easy 
in this case but in general can be quite tedious (though completely algorithmic) . 
The second is finding the solution to a linear system, which was used to solve 
for n and m in the G step, and used to solve for k in the G2 step (in a free 
group with more generators, this would have been a linear system with more 
than one equation). 

As is to be expected, computation of twisted conjugacy classes in many cases 
may require checks in G n for n > 2. Such examples can easily be constructed 
by a computer search. 

Example 4.3. Let G — (a,b), and let 

a 1— » aba -1 

V- b „ a -2 b 4 

Theorem 12.11 gives 

RT{tp) = piaba- 1 - a~ 2 - a~ 2 b - a~ 2 b 2 - a~ 2 b 3 ). 

It can be verified that all five of the above terms are twisted conjugate in G n for 
n G {1,2,3}, but none are twisted conjugate in G4. Thus the Nielsen number 
is 5. 
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We now turn to the question of verifying equality of twisted conjugacy 
classes. Applying the process described above to two words which are in fact 
twisted conjugate will result in a non-terminating sequence of computations in 
the groups G n , each time resulting in solutions for the various elements above 
labeled z. 

To avoid such an infinite computation, we propose a technique in the spirit 
of Remark 14.21 In the abelianization G, if a solution element z is obtained, a 
sequence of "candidates" for testing twisted conjugacy in G is constructed by 
producing all possible reorderings of the generators appearing in z. Thus if we 
obtain an element z — a — 2b, our list of candidates will be 

a&- 2 ,& _1 a& _1 ,6- 2 a. 

Each of these elements can be tested by twisted conjugation as a candidate for 
realizing the twisted conjugacy in G. 

A similar process can be carried out in G n for n > 1: after obtaining z S G n , 
our list of candidates is obtained by inserting the weight n basic commutators 
appearing in z in all possible orderings and in various forms into each of our 
candidates previously obtained in our check from G n _i. 

Example 4.4. Let G = (a,b), and 

a i — > a 2 ba 
V b h-v b 2 a 

We will use the above candidates checking construction to show that [a] = [a 2 b] 
(the elements a 2 b and a appear in RT(<p)). 

Our check in G shows that the two elements are twisted conjugate by —b. 
Thus our only candidate from G is the element but a computation shows 
that 

y>(6 -1 )(o 2 6)6 = a^ZrVfc 2 ^ a. 

Now we do a check in G%, and we find that the two elements are twisted conju- 
gate by & _1 [a, b]^ 1 E Gi- Now our list of candidates is: 

,9i = 1 [a, £>] 1 = ab- x ar x 

g 2 = fc-^crSft -1 ] -1 = fc- 2 a _1 6a 

g 3 = 6 -1 [cr\&] = b^a^bab- 1 

g A = b-ifab- 1 ] = b- 1 ab- 1 a~ 1 b 

g 5 = [a,6] -1 6 _1 = bab^a^b' 1 

g 6 = [a-^b- 1 ]-^- 1 = 6- 1 o~ 1 6o6 -1 

g 7 = [o _1 ,6]6 _1 = orHab' 2 

g$ = [a, = a6 _1 a _1 

Computing each of ( <7i )a 2 6(<7i) — 1 reveals that ip(gi)a 2 bgi 1 — a, and so [a] = 
[a 2 b]. 
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Note that in the above example we inserted the basic commutator [a, 6] _1 
into the word b^ 1 in each of two positions in one of four forms, these being the 
four versions of this commutator in G which become [a, when projected into 
G2. Constructing the various forms of a weight 3 commutator to use in a list of 
candidates would be cumbersome, and we do not attempt this construction for 
nilpotency class higher than 2. 

As an alternative to the candidates construction procedure, a more pedes- 
trian approach is always available: after a check for twisted conjugacy in G n , 
use as list of candidates all words in G of word length n. This produces a dif- 
ferent list of candidates from that described above, but is guaranteed to find an 
element realizing the twisted conjugacy if n is sufficiently high. 

5 Success rates 

The technique exhibited in the examples above can be summarized as follows: 
Starting with n = 1, write an expression for a generic element z G G n and 
the element (p(z) in Hall normal form. Solve a linear system to decide if the 
elements are twisted conjugate in G n . If the elements are not twisted conjugate 
in G ni then the elements are not twisted conjugate in G. If the elements are 
twisted conjugate in G n by a unique element z, then construct a finite list of 
candidates (either intelligently by using the structure of z or in a brute force 
manner by taking all words of length n) for testing twisted conjugacy in G. If 
all of these candidates fail, then increment n and repeat the above. 

This technique will decide any given twisted conjugacy problem provided 
that the following statement is true: If cp : G — > G is a map on a free group, 
and g and h are two elements of G which are not twisted conjugate in G, then 
there is some n for which g and h are not twisted conjugate in G n . Thus any 
non-twisted-conjugate elements will be detected as such in G n for some n. Such 
a statement does not hold in general, though, as the following argument shows. 

Proposition 5.1. For any <p : G — > G, if g, h E G are words such that (p n (g) G 
h-f n (G) for all n, then [g] = [h] in G n for all n. 

Proof. The proof is based on the fact that [y(x)] = [x] for any element x, since 
ip(x) — ip(x)xx~ 1 . Iteration of the map gives [<p"(a;)] = [x] for any n. 

Now for our element g, we have ^5™(g) = h in G n , and so in particular 
[£>"(?)] = [h]- But = [g] by the above, and so we have [g] = [h]. □ 

We can use the above to build homomorphisms ip : G —>jG with the property 
that there are words g,h G G with [g] ^ [h] but [g] = [h] in every nilpotent 
quotient G n . 

Example 5.2. Let G — (a, b), and let (p : G — > G be the map 

a 1 — ► [6, a] 
b » a- x b 
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Now we have ip(j n (G)) C j n +i(G) for n > 0, since <p replaces a with a weight 
2 commutator and all commutators in G involve the element a. Since <p(a) G 
71(G), we have <£ n (a) S 7«(G) for all n, and thus that [a] = [1] in G n for all n. 

But <y5 is a mapping with remnant, and Wagner's algorithm can be used to 
show that in fact [a] ^ [1] in G. 

Though the above example shows that nilpotent quotients cannot be used 
directly to solve twisted conjugacy problems in all cases, the technique gives very 
good success rates in experimental testing. We have created an implementation^ 
of the process in the computational algebra system Magma | Bos ma et al, 1994] . 

There are two ways that an implementation of this tcchinque can fail to 
decide any given twisted conjugacy problem. One type of failure is when the 
linear system which arises in the twisted conjugacy computation in G n has 
infinitely many solutions. In such a case the check in G n +i will require finding 
an integral solution of a polynomial system, which will in general be difficult. 
Such a failure can only occur when the coefficients in the linear system give a 
singular matrix, which we expect to occur relatively infrequently. 

A second type of failure is that the implementation exhausts its resources 
in computing the required Hall normal forms in G n . Even for groups of 2 or 3 
generators, human computation of the Hall form in class 3 or 4 is barely feasible. 
Since the number of basic commutators of weight n grows exponentially in n 
we expect that any computer implementation could conceivably exhaust its 
resources before detecting that two given twisted conjugacy classes are indeed 
distinct. We expect this type of failure to occur increasingly in free groups with 
large numbers of generators. 

Tables [1] and [2] give experimental results of application of the nilpotent quo- 
tients technique to 10,000 randomly generated mappings on the free group on 
k generators for k = 2,3,4. In each number I 6 {2,3,4,5} is chosen 

and a mapping is generated by assigning the image of each generator to be a 
randomly chosen word of length at most I. Theorem 12. II is applied to give a list 
of group elements which must be divided into their twisted conjugacy classes. 
A sucessful computation is one which is able to decide the twisted conjugacy of 
these elements. 

Table [T] gives the rates of each of the two types of failure above. A "ma- 
trix failure" is declared when the linear system computation results in infinitely 
many solutions. A "complexity failure" is declared when the nilpotency class 
reaches 5, as this is the level at which the computation of the Hall normal form 
becomes difficult for our implementation (run on a personal computer, current 
in 2005). Since a single mapping can trigger both types of error (computa- 
tion of RT(ip) in general requires several twisted conjugacy decisions), the row 
percentages may not sum to 100%. The column labeled "average depth" gives 
the average nilpotency class required to distinguish twisted conjugacy classes 
in these random mappings. A depth of n indicates that a check in G n was 
necessary. The column labeled a gives the standard deviations of the depths. 

1 Available from the author's web site: http://www.messiali.edu/~cstaecker 
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k 


I 


Success 


Matrix failure 


Complexity failure 


Avg. depth 


a 




2 


94.27% 


4.30% 


2.97% 


1.10 


0.28 




3 


90.46% 


6.66% 


4.09% 


1.13 


0.34 


2 


4 


85.91% 


7.75% 


9.06% 


1.18 


0.42 




5 


84.20% 


9.89% 


8.66% 


1.17 


0.38 




2 


87.45% 


11.26% 


4.38% 


1.18 


0.45 


3 


3 


85.60% 


13.77% 


4.74% 


1.17 


0.42 


4 


85.20% 


12.38% 


6.19% 


1.21 


0.47 




5 


84.46% 


13.36% 


6.29% 


1.23 


0.45 




2 


88.87% 


10.95% 


3.48% 


1.13 


0.41 




3 


85.60% 


13.77% 


4.74% 


1.16 


0.42 


4 


4 


84.51% 


14.56% 


5.87% 


1.21 


0.45 




5 


84.13% 


15.04% 


5.64% 


1.21 


0.42 



Table 1: Results of testing for success rates on random mappings of word length 
I on the free group on k generators. 



n 


I 


Nilpotent quotients 


Abelianization 


Wagner's Alg. 




2 


94.27% 


82.75% 


41.80% 




3 


90.46% 


70.45% 


48.32% 


2 


4 


85.91% 


59.07% 


54.71% 




5 


84.20% 


51.14% 


59.83% 




2 


90.98% 


77.96% 


15.54% 




3 


87.45% 


65.99% 


24.70% 


3 


4 


85.20% 


55.17% 


33.72% 




5 


84.46% 


47.43% 


41.45% 




2 


88.87% 


76.50% 


5.84% 




3 


85.60% 


64.76% 


13.67% 


4 


4 


84.51% 


54.63% 


22.63% 




5 


84.13% 


47.81% 


30.18% 



Table 2: Comparison of success rates of various techniques on random mappings 
of word length I on the free group on k generators. 
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Table[2]gives our sucess rates compared to the existing techniques of abelian- 
ization and Wagner's algorithm. The technique used for data in the "Abelian- 
ization" column uses only the abelianization for distinguishing classes, and uses 
only the identity [f(g)] = [g] for equating classes (this is the general strategy 
employed in |Hart, 20051 ). The column labeled "Wagner's alg." records the 
percentages of maps satisfying Wagner's remnant condition, for which her algo- 
rithm will apply (that the percentages grow in I is expected in light of Theorem 
3.7 of | Wagner, 1999p . 



6 Surfaces without boundary 

The nilpotent quotients process can be used with minimal modifications when 
G is the fundamental group of a compact hyperbolic surface without boundary. 
Work of Fadell and Husseini in |Fadell and Husseini, 1983| together with a tech- 



nique by Davey, Hart, and Trapp |Davey et ai, 199"6 reduce the computation 



of the Nielsen number on compact hyperbolic surfaces without boundary to the 
computation of twisted conjugacy classes in the fundamental group. 

Wagner's technique does not apply if G is not a free group, and neither will 
the techniques of [Bogopols ki et ai, 2006| . Because surface groups are easily 
expressed in terms of commutator relations, we can use the nilpotent quotients 
technique with only trivial modifications in this setting. 

The only modification that must be made is to the precise structure of the 
Hall normal form. For instance, if G is the fundamental group of the genus 2 
compact surface, then G has group presentation G = (a, b, c, d\ [a, b] [c, d\ = 1). 
The Hall normal form of an element of (e.g.) G^ can be obtained by applying 
commutation rules as if G were the free group on 4 generators, along with an 
additional rule that [c, d]^ 1 = [a,b]. It is convenient for us that the group 
structure of G is so compatible with the Hall normal form. 

The case of surfaces without boundary is somewhat more difficult to imple- 
ment in Magma, as it involves computations in finitely-presented rather than 
free groups. The capabilities of Magma are somewhat lacking in this regard- 
in particular Magma (as of version 2.13-15) is unable to reliably solve the word 
problem in a surface group (although this word problem is solvable). This causes 
the candidates checking process to return false negatives, as the implementa- 
tion may not recognize when two elements are actually equal. Statistics such as 
those in Table Q] are also difficult to produce in this setting as it is difficult to 
generate random endomorphisms of surface groups. 



7 Doubly twisted conjugacy 

We conclude with a brief discussion of how our technique can be applied to 
the doubly twisted conjugacy relation: Given two maps ip,ip : G — > H and two 
elements h, k G H, we say that h and k are (doubly) twisted conjugate (we write 
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[h] = [k]) if there is an clement g G G with 

h = ip{g)ktp(gy 1 . 

This relation is fundamental in Nielsen coincidence theory (see |Gongalves, 2005| ), 
playing the same role as ordinary twisted conjugacy in fixed point theory. 

For any n, the maps cp and ip will induce maps ip,ip : G n — > H n , and 
the doubly twisted conjugacy relation can in principle be solved by using Hall 
normal forms in G n and H n just as in the ordinary twisted conjugacy problem. 

Example 7.1. Let G = H = (a,b), and let our maps be 

a i— > b 2 a , a i— > a 3 
^ : 6 ^ a- 2 ^ : 6 ^ a- 1 

We will decide the twisted conjugacy of the elements b and 6 _1 . 

We begin with check in the abelianization, where any element z € G has the 
form z = na + mb. We compute that <p(z) = (n — 2m)a + 2nb and —^(z) = 
(— 3n + m)a, and thus we have 

tp(z) - b- tfj(z) = (— 2n - m)a + (2n - 1)6. 

Equating this with 6 and solving gives n = 1 and m = — 2. This solution in the 
abelianization gives three candidates for twisted conjugacy: 

a6- 2 ,&- 1 a&- 1 ,6- 2 a, 

but checking each shows that none of these realize the twisted conjugacy in H. 

We proceed to the class 2 nilpotent quotient, where any element z £ Gi has 
the form z = a n b m \a, b] k . Our computation in G shows that n = 1 and m = —2, 
simplifying our element to z = ab~ 2 \a, 6] . We compute 

£(z) =6 2 S(S- 2 )- 2 [6 2 S,a- 2 ] fe = a 5 6 2 [S 5 ,6 2 ][6 2 ,a- 2 ] fc =a 5 b 2 [a,b] w+ik , 
4>(z) =a 3 (a- 1 )- 2 [a 3 ,a- 1 ] k =a 5 , 

and so 

^Jft- 1 ^)- 1 =a 5 b 2 [a,b] w+ik b- 1 a- 5 = a 5 ba- 5 [a,b] w+ik 
= fe[a- 5 ,&][a,6] 10+4fc - fe[a,6] 5+4fe . 

Equating this with b gives 5 + 4fc = which is impossible for integral k. Thus 

[fc^pr 1 ]. 

There is in the literature no analogue of Theorem 12.11 in coincidence theory, 
but presumably one may be available in the future, and our technique is cur- 
rently the only available technique for distinguishing doubly twisted conjugacy 
classes (no version of Wagner's algorithm is known in coincidence theory, and 
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32.88% 
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98.50% 


1.15% 
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Table 3: Success rates for doubly twisted conjugacy relations. Random map- 
pings of word length 3 from the free group on k\ generators to the free group on 
ki generators were tested in deciding twisted conjugacy between two random 
words of length at most 3. 



the methods of Bogopolski et ai, 2006| do not extend in an obvious way to 
doubly twisted conjugacy). 

Table [3] gives success rates for the technique applied to 10,000 randomly 
generated twisted conjugacy relations. In each case, two random mappings of 
"word length" 3 (the quantity labeled I in Tables Q] and \2§ are generated from 
the free group on k\ generators to the free group on ki generators. Two random 
elements of the codomain group are generated with word length at most 3, and 
the implementation attempts to determine their twisted conjugacy. Entries in 
the table with no digits to the right of the decimal point are exact figures, e.g. 
in the case where k\ = 4 and k^ = 3 the depth was exactly 1 in each of the 
10,000 test cases, and there were exactly complexity failures. 

Note that the technique is much less successful if the rank of the domain is 
greater than the rank of the codomain. This is to be expected, as the technique 
will fail when our linear system computation (always having more variables than 
equations if fc x > fc 2 ) yields infinitely many solutions. Note that such cases are 
not handled by our Magma implementation, but could in principle be done by 
hand. These would require finding integer solutions to polynomial systems, and 
so we do not always expect the computation to be successful, but particular 
examples may be computable. 

Especially striking are the extremely high sucess rates when > fci. The 
vast majority of these twisted conjugacy relations are decided in the abelianiza- 
tion (and with negative result), as the overdetermined linear systems are unlikely 
to have any solutions. For example in the case of k\ =3 and &2 = 5, an equiva- 
lence in the abelianization would require our random elements to satisfy a linear 
system of 5 equations and 3 unknowns. If the elements are indeed equivalent 
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in the abelianization, a further equivalence in the class two nilpotent quotient 
would require our random data to satisfy a linear system of 10 equations and 3 
unknowns (the free group on 5 generators has 10 basic weight 2 commutators, 
and the free group on 3 generators has 3 basic weight 2 commutators). This 
is not, of course, to say that twisted conjugate elements do not occur in such 
cases where &2 > k\, but only that this occurs very infrequently in the case of 
two words of length 3. 

We have omitted from our testing the cases where k\ or fc 2 is 1. The twisted 
conjugacy relation is generally solvable in these cases by hand. 

Let ip,tp : G —* H are maps of free groups where G — (a), and let h,k £ H 
be two words. Then the twisted conjugacy problem is equivalent to finding some 
integer n with 

1 = /i-V(a)"fcV(a)" n , 

and this can typically be confirmed or denied by inspection. 

If tp, ip : G — ► H are maps of free groups and H has rank 1, let G = 
(oi, . . . ,a n ), and note that <^(7i(G)) C 71 (if) = 1 since H is abelian. Thus 
we have ip(uv) — <p(vu[v,u]) = ip(vu) for any words v,u £ G, and similarly 
tp(uv) = tp(vu). Thus, though there is no convenient normal form for elements 
z £ G, we can say in general that 

<p{z) = p(ar ■ ■ ■ <C") ; 1>{z) = ^ ■ ■ ■ m*») 

by rearranging the generators of z. 

Now to decide the twisted conjugacy of elements h,k £ H, we examine the 
equation 

h = <p(aF...<C~)ty(ai m K..a; m ~), 

and we will be able to determine whether or not the above has solutions because 
it is an equation in the abelian group H. If the above has a solution, then the 
elements are twisted conjugate, and if not, they are not. 
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